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1  SUMMARY 


The  research  objective,  to  develop  new  objective  tools  to  aid  in  software  verification  and  validation 
based  on  topological  methods,  found  that  the  application  of  topology  to  software  with  composite 
non-numerical,  non-metric  data-types,  proved  to  be  too  course-grained  to  provide  a  mathemati¬ 
cal  foundation  for  rigorously  reasoning  about  the  evolution  of  the  composite  state  of  executing 
software  but  overcame  this  difficulty  by  applying  a  refinement  of  topological  structure,  called  con¬ 
vergence  structure  (CS). 

CS  was  applied  to  serve  as  the  mathematical  underpinning  for  specifying  the  evolution  of 
composite  states  by  providing  an  analytical  foundation  for  the  evolution  of  such  states.  Such 
a  foundation  is  geometric  in  character  and  is  derivable  from  symmetries  in  the  state  space  by 
group-theoretic  means  and  then  determining  optimal  CS  with  respect  to  which  the  symmetries  are 
continuous.  This  approach  is  justified  because  the  evolving  values  of  variables  that  take  entire 
data- structures  as  their  values  cannot  in  general  be  numerically  encoded  without  either  having 
to  subjectively  assess  the  appropriateness  of  a  particular  encoding,  or  else  encounter  question¬ 
begging  circularities.  At  the  logical  level  in  an  actual  computing  device,  composite  datatype  values 
are  networks  strewn  through  the  device’s  memory;  think  of  a  variable  whose  type  is  an  array  of 
balanced  trees.  Some  means,  structural  in  nature,  is  required  to  rigorously  specify  the  evolution  of 
composite  states  involving  non-numerical,  non-metric  components.  Composite  datatypes  cannot 
in  general  be  numerically  encoded  without  either  having  to  subjectively  assess  the  appropriateness 
of  a  particular  encoding,  or  else  encounter  question-begging  circularities.  Since  the  fundamental 
aim  of  this  project  has  been  to  objectify  correctness  assessments  of  software  behavior,  the  project 
required,  at  the  very  least,  a  means  for  specifying  correct  behavior  of  the  evolution  of  a  composite 
state  having  non-numerical  components. 

The  verdict  as  to  whether  the  application  of  persistent  homology  to  the  problem  of  obtaining 
objective  signatures  that  would  indicate  the  present  or  absence  of  critical  behaviors  of  execut¬ 
ing  software  remains  undecided  at  the  end  of  the  project  period,  while  unsupported  work  on  this 
problem  continues.  During  the  project,  in  cooperation  with  Dr.  W.  A.  Miller  of  Florida  Atlantic 
University  (FAU),  Dr.  Matthew  Come,  (NRC  Fellow  at  AFRL/RI),  doctoral  student  Shannon  Ray 
(FAU),  and  others,  a  computational  experiment  was  conducted  consisting  of  the  application  of 
persistent  homology  to  a  2-dimensional  rotational  solid  that  collapses  and  to  3 -dimensional 
dumbbells  that  manifest  neckpinch  singularities.  The  appearance  of  critical  geometric  phenom¬ 
ena  in  these  models  was  compared  with  the  results  produced  by  persistent  homology  methods  and 
we  concluded  via  various  resolutions  that  persistent  homology  does  indicate  geometric  criticality. 
However,  these  conclusions  were  reached  by  ad  hoc  means  and  an  objective  methodology  to  reli¬ 
ably  determine  even  the  general  kinds  of  dynamical  geometric  phenomena  indicated  by  persistent 
homology  remains  to  be  developed.  Moreover,  the  geometric  phenomena  associated  with  neck¬ 
pinch  singularities  occurring  in  a  collapsing  dumbbell  is  vastly  simpler  than  the  flow  through  a 
state  space  produced  by  executing  software.  Moreover  in  the  course  of  the  project  an  additional 
computational  experiment  was  performed  with  the  assistance  of  consisting  of  applying  persistent 
homology  to  detect  changes  to  the  structure  of  a  small  but  nontrivial  Petri  net  for  the  individual 
steps  of  any  given  trajectory  could  have  been  calculated  by  hand,  but  for  which  the  persistent 
homology  calculations  required  a  computer.  This  experiment  was  carried  out  with  the  assistance 
of  Dr.  Konstantin  Mischaikow  (Rutgers  University),  Dr.  Vidit  Nanda  (then  of  Rutgers 
Univer-sity),  Mr.  Anthony  Ford  (then  of  AFRL/RI)  and  Dr.  Jonathan  McDonald  (then  an  NRC 
Fellow 
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at  AFRL/RI).  The  ad  hoc  result  was  that  gross  changes  in  the  dynamics  of  the  Petri  net  under 
various  loads  could  be  detected,  but  more  subtle  changes  in  the  dynamics  could  not.  This  does 
not  deprecate  persistent  homology  as  an  approach;  rather  it  indicates  that  a  vast  amount  of  close 
investigation  of  the  relationship  between  the  Betti  spectra  as  a  function  of  a  filtration  parameter 
produced  by  persistent  homology,  and  the  dynamically  evolving  underlying  geometry  is  needed, 
and  will  best  be  served  by  understanding  much  more  about  the  dynamics  of  such  flows,  first. 

Accordingly,  the  project  focused  on  developing  the  mathematically  rigorous  specification  and 
verification  of  dynamical  systems  involving  variables  with  compound  type.  We  call  these  dynam¬ 
ical  systems  heterotic  dynamical  systems.  It  is  easy  to  see  a  program  with  variables  of  composite 
type  as  a  heterotic  dynamical  system,  and  it  is  to  be  expected,  therefore,  that  there  will  be  com¬ 
ponents  of  the  evolving  computation  states  associated  with  the  execution  of  these  programs  whose 
type  is  essentially  a  class  of  discrete  structures  such  as  linked  lists  for  example. 

The  capability  to  rigorously  define  heterotic  dynamical  systems  and  consequently  a  specifica¬ 
tion  logic  for  them  is  due  to  key  properties  of  the  collection  of  all  convergence  spaces.  Convergence 
spaces  and  the  continuous  functions  among  them  form  a  large  Cartesian-closed  Category ,  which 
is  simultaneously  a  Cartesian  closed  completion  of  the  category  of  topological  spaces  and  the  cat¬ 
egory  of  reflexive  directed  graphs  whose  arrows  are  directed  graph  homomorphisms.  Within  this 
category  the  conditions  for  the  functions  to  among  the  convergence  spaces  are  the  same:  preser¬ 
vation  of  converging  filters.  Thus,  homorphisms  and  continuous  functions  are  revealed  to  be  the 
same  kind  of  mathematical  object:  a  function  that  preserves  convergence  of  filters. 

A  doctoral  student,  D.  R.  Patten,  investigated  CS  methods  on  discrete  structures  for  the  pur¬ 
pose  of  rigorously  formulating  heterotic  dynamical  systems  involving  discrete  structure  data-types. 
Patten’s  investigation  will  aid  in  developing  formal  specifications  for  object-oriented  and  similar 
software  involving  composite  data-types.  Patten  successfully  defended  his  doctoral  dissertation 
embodying  this  research  in  June,  2014.  Patten’s  work  identified  a  substantial  amount  of  practical 
detail  to  exploit  when  constructing  differential  calculus  among  spaces  with  components  that  are 
discrete  structures  such  as  networks. 

2  INTRODUCTION 

The  application  of  topology  to  verifying  and  validating  software  with  composite  non-numerical, 
non-metric  data-types,  is  too  course-grained  to  provide  a  mathematical  foundation  for  rigorously 
reasoning  about  the  evolution  of  the  composite  state  of  executing  software.  A  topology  on  a  set  of 
elements  is  about  nearness:  which  elements  are  near  to  which  elements.  The  individual  elements 
can  have  a  complicated  internal  structure,  and  there  may  not  be  enough  separation  strength  [20] 
to  numerically  characterize  the  nearness  notion  with  a  metric  or  with  continuous  maps  into  the 
real  numbers  [42].  Accordingly,  applications  of  computational  algebraic  topology  that  face  such 
limitations  require  appropriate  adaptations  of  topological  methods.  The  objective  of  this  research 
therefore  became  that  of  adapting  techniques  of  computational  algebraic  topology  to  much  finer- 
grained  topology-like  convergence  structures  of  which  topologies  are  a  special  case  [10]. 

We  know  that  such  an  adaptation  can  be  carried  out  for  the  following  reason:  A  fundamental 
motive  for  topology  is  to  capture  the  structural  essentials  of  the  notion  of  continuity  shorn  of  inci¬ 
dental  features,  such  as  quantitative  characterizations  of  continuity  familiar  from  the  e,  5  definitions 
of  elementary  calculus.  Nevertheless,  the  relatively  coarse-grained  nature  of  topologies  seems  to 
be  unsuitable  for  the  analogous  characterization  of  differentiation,  while  convergence  structures  do 
allow  for  a  structural  non-quantitative  characterization  of  differentiation  seamlessly  with  classical 
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analysis  [14].  There  are  several  reasons  underlying  this  capability  of  convergence  spaces',  i.e.  a 
convergence  space  is  a  nonempty  set  equipped  with  a  convergence  structure. 

1.  The  fine-grained  character  of  convergence  structures:  The  distinction  between  convergence 
structure  and  topological  structure  can  be  succinctly  grasped  by  considering  the  specializa¬ 
tion  of  both  to  directed  graphs  with  self-loops,  i,e,  with  a  reflexive  edge  relation.  A  directed 
graph  with  a  reflexive,  and  transitive  edge  relation  is  a  topological  space  where  every  vertex 
occurs  in  a  smallest  open  set  containing  that  vertex  and  consisting  of  all  adjacent  vertices.  A 
directed  graph  with  a  merely  reflexive  edge  relation  is  a  convergence  space.  Thus  by  coars¬ 
ening  a  directed  graph  by  transitively  closing  its  edge  relation  one  moves  from  a  convergence 
structure  to  the  induced  topological  structure. 

2.  The  Cartesian-closure  of  the  category  of  convergence  spaces:  In  a  Cartesian  closed  category 
[33],  among  other  features,  each  collection  of  arrows  Hom(A,  B)  from  an  object  A  to  an  ob¬ 
ject  B  is  represented  in  the  category  by  at  least  one  object.  The  collection  of  all  convergence 
spaces  forms  a  Cartesian-closed  category  CONV  in  which  the  objects  are  convergence  spaces 
and  the  arrows  are  continuous  functions  among  them  [24]  [17].  As  objects  of  the  category, 
the  function  spaces  Yx  of  continuous  functions  from  X  to  Y  are  themselves  convergence 
spaces  with  a  convergence  structure  that  is  determined  by  the  convergence  structures  on  X 
and  Y  in  a  uniform  manner  throughout  CONV.  CONV  is  a  cartesian-closed  extension  of 
both  the  category  RDIGRAPH,  the  category  of  directed  graphs  with  self-loops  where  di¬ 
rected  graph  homomorphisms  are  the  arrows,  as  well  as  the  category  TOP  of  topological 
spaces  with  continuous  functions  as  the  arrows.  RDIGRAPH  and  TOP  axe,  full  subcategories 
of  CONV,  i.e.  no  arrows  are  lost  in  restricting  to  these  subcategories. 

3.  A  differential  calculus  on  a  collection  of  spaces:  In  prior  work  the  PI  formulated  differential 
calculus  on  convergence  spaces.  [14].  A  differential  calculus  is  determined  by  the  relation¬ 
ship  that  function  g  is  a  differential  of  function  /  at  a  point  x,  together  with  a  functorial  chain 
rule.  In  other  words,  the  (covariant)  derivative  operator  maps  points  in  the  domain  of  /  to 
differentials.  The  functoriality  of  the  chain  rule  is  what  yields  calculus;  for  example  the  other 
familiar  rules  of  elementary  calculus  on  Euclidean  spades,  such  as  the  product  rule  follow 
from  the  chain  rule  together  with  differentials  of  sum,  product  and  projection,  and  the  rules 
for  cartesian  products  of  functions.  To  obtain  a  differential  calculus  using  non-quantitative, 
non-metric  structural  means  on  collection  of  spaces,  a  functorial  chain  rule  must  be  obtained, 
and  functions  among  convergence  spaces  must  be  selected  to  serve  as  differentials.  For  a  cal¬ 
culus  to  be  usefully  robust,  we  prefer,  but  do  not  require,  that  the  projection  functions  and 
diagonal  functions  such  as 


(x,  y)  i-»  x  and  x  i-»  (x,  x) 

have  suitable  differentials.  The  continuity  of  these  sorts  of  functions  having  to  do  with 
Cartesian  products,  as  well  as  continuity  of  the  composition  operator  that  allows  to  combine 
two  functions  in  sequence,  are  the  features  of  the  Cartesian-closedness  of  CONV  essential 
to  our  purpose. 

4.  The  focus  on  computational  algebraic  topology:  To  date  this  is  largely  based  on  homology 
for  reasons  of  computational  tractability.  The  simplex  is  the  basic  mathematical  object  at  the 


Approved  for  Public  Release;  Distribution  Unlimited. 

3 


core  of  mainstream  homology.  The  standard  geometrical  representation  of  an  n-dimensional 
simplex  is  the  convex  hull  of  n  + 1  orthonormal  basis  vectors  in  n  + 1  dimensional  Euclidean 
space.  However,  in  a  more  general  topological  space  X,  even  a  metric  space,  a  simplex  is  a 
continuous  mapping  from  a  standard  geometrical  representation  simplex  into  X.  The  target 
space  of  such  a  continuous  map  can  just  as  well  be  a  convergence  space,  i.e.  a  nonempty 
set  equipped  with  a  convergence  structure,  and  must  be,  if  the  target  space  is  a  dynamical 
system’s  state  space  with  components  of  the  states  that  consist  of  discrete  data-structures. 
Since  RDIGRAPH  occurs  in  multiple  guises  as  a  full  subcategory  of  CONV  in  a  special  way: 
there  are  continuous  functions  among  reflexive  digraphs;  i.e.  there  are  continuous  functions 
among  discrete  structures.  This  shows  that  homomorphisms ,  are  a  special  case  of  continuous 
functions.  This  observation  does  not  hold  within  the  category  of  topological  spaces,  and  this 
is  where  CONV  shines:  The  set  YY  of  all  continuous  functions  on  Y  is  also  a  convergence 
space,  where  if  Y  is  a  reflexive  digraph,  then  Yy  is  the  set  of  all  endomorphisms  on  Y . 
Y]  has  subspaces  consisting  of  abelian  groups  of  automorphisms  of  Y .  These  subspaces 
have  group  endomorphisms  on  them  from  which  arise  ring  modules  that  are  semi-regular 
actions  on  Y,  decomposing  Y  into  orbits  of  the  action.  The  endomorphisms  in  the  module 
are  the  translations  of  Y  in  exactly  the  same  sense  as  affine  functions  on  Euclidean  spaces 
are  translations  of  those  spaces.  Since  these  translations  are  also  automorphsisms  of  Y  they 
are  symmetries  of  Y .  These  ring  modules  are,  in  effect,  the  natural  geometries  inherent  in  the 
convergence  structure  on  Y.  Whenever  the  ring  modules  are  algebraic  fields,  the  individual 
orbits  into  which  Y  decomposes  are  vector  spaces.  Thus  Y  is  revealed  as  a  disjoint  union 
of  vector  spaces  in  such  cases,  but  is  always  a  disjoint  union  of  ring  modules  in  any  case. 
Ring  modules  may  be  only  nearly  vector  spaces  but  for  that  fact  that  not  all  of  the  nonzero 
scalars  need  to  be  invertible  and  arise  naturally  in  functional  analysis  as,  for  example,  the 
space  of  differential  forms  on  a  smooth  manifold.  Vector  spaces  are  special  cases  of  ring 
modules.  The  work  in  this  project  culminated  in  the  study  of  the  persistent  homology  to  a 
2-dimensional  rotational  solid  that  collapses  and  to  3-dimensional  dumbbells  that  manifest 
neckpinch  singularities.  This  study  indicated  that  applications  of  persistent  homology  to 
high  dimensional  complicated  dynamical  systems  is  premature  [5,  15]. 

The  project  focused  on  developing  the  mathematically  rigorous  specification  and  verification  of 
dynamical  systems  involving  variables  with  compound  type.  We  call  these  dynamical  systems 
heterotic  dynamical  systems.  There  is  a  very  large  degree  of  freedom  in  the  choice  of  conver¬ 
gence  structure  with  which  to  equip  a  state  space.  This  range  of  choice  must  be  severely  narrowed 
if  pre-existing  software  is  to  be  modeled  by  heterotic  dynamical  systems  in  an  objective  sense. 
Convergence  structure,  like  topological  structure  but  more  sensitively,  determines  in  a  structural 
sense,  what  elements  that  are  near  to  each  other.  Specific  executions  of  software  produce  trajec¬ 
tories  in  their  computation  space  that  can  be  regarded  as  continuous  images  (remember  that  with 
convergence  spaces  our  naive  intuitions  about  what  continuity  is  often  fail,  and  this  is  true  even 
with  conventional  applications  of  topology)  of  time,  either  modeled  as  the  integers  in  the  case  of 
fully  autonomous  digital  computation,  or  as  the  real  numbers  if  the  software,  for  example,  is  part 
of  a  control  system  with  feedback.  The  execution  trajectories  determine  a  finest  (and,  therefore, 
contravariantly  a  smallest)  convergence  structure  on  the  state  space.  The  point  is  that  the  execution 
trajectories  uniquely  determine  the  CS,  and  hence  the  geometric  structure  on  the  state  space,  and 
conversely,  for  a  heterotic  dynamical  system  defined  relative  to  an  absolute  time,  the  convergence 
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structure  on  the  state  space  constrains  the  possible  trajectories. 

A  doctoral  student,  D.  R.  Patten,  [36]  investigated  CS  methods  on  discrete  structures  for  the 
purpose  of  rigorously  formulating  heterotic  dynamical  systems  involving  discrete  structure  data¬ 
types.  Patten’s  investigation  will  aid  in  developing  formal  specifications  for  object-oriented  and 
similar  software  involving  composite  data-types.  Patten  successfully  defended  his  doctoral  dis¬ 
sertation  embodying  this  research  in  June,  2014.  Patten’s  work  identified  a  substantial  amount  of 
practical  detail  to  exploit  when  constructing  differential  calculus  among  spaces  with  components 
that  are  discrete  structures  such  as  networks. 

3  METHODS,  ASSUMPTIONS  AND  PROCEDURES 

3.1  Differential  calculus  on  nonstandard  spaces 

With  topology,  continuity  of  functions  generalizes  from  the  context  of  classical  analysis  to  a  huge 
collection  of  structures,  the  topological  spaces,  ity.  Differentiability  also  generalizes  to  all  topolog¬ 
ical  spaces,  and  beyond  (as  does  continuity)  to  all  convergence  spaces.  Just  as  continuity  neither 
presupposes  any  seperation  strength  nor  any  notions  of  linearity,  neither  does  differentiability.  The 
familiar  differential  calculus  on  Euclidean  spaces  is  of  course  intrinsically  dependent  on  the  vector 
space  structure,  but  this  is  due  to  the  choice  of  functions  used  to  serve  as  differentials,  and  the 
consequent  determination  of  the  conditions  under  which  functions  are  differentiable.  What  counts 
is  the  relation  “differential  g  is  a  differential  of  /  at  x”.  Unless  we  demand  of  g  that  it  satisfy 
some  kind  of  linearity  property,  linearity  does  not  intrinsically  enter  into  the  relation.  CONV,  the 
Cartesian-closed  category  of  convergence  spaces  and  continuous  maps  [30],  is  central  in  this  gen¬ 
eralization.  The  spaces  and  functions  of  interest  are  naturally  organized  into  categories  and  it  is 
important  to  note  the  nature  of  the  containments  and  embeddings  that  are  involved.  For  example, 
any  method  for  constructing  a  differential  calculus  for  mappings  between  arbitrary  convergence 
spaces  gives  such  a  method  for  all  reflexive  digraphs  and  all  topological  spaces. 

The  category  CONV  includes  all  topological  spaces  but  enjoys  several  substantial  advantages 
over  the  category  TOP  of  topological  spaces.  Importantly  for  computation,  CONV  contains  all 
reflexive  directed  graphs  (finite  and  infinite).  Among  digraphs,  continuity  is  the  property  of  be¬ 
ing  edge-presen’ing,  i.e.  a  digraph  homomorphism.  But,  powerfully,  and  unlike  TOP,  CONV  is 
a  Cartesian-closed  category.  Several  immediate  consequences  of  Cartesian-closure  and  the  rela¬ 
tionship  between  TOP  and  CONV  are:  (1)  convergence  spaces  preserve  the  notion  of  continuity 
on  topological  spaces;  (2)  convergence  spaces  allow  fine  control  over  continuity,  and  in  various 
circumstances  permit  us  to  strengthen  the  conditions  for  continuity;  (3)  there  is  a  uniform  way 
of  regarding  all  spaces  of  continuous  functions  as  convergence  spaces  (enabling  hybrid  functional 
analysis)  and  what  is  key  to  differential  calculi:  (4)  the  composition  operation  on  finite  products  of 
function  spaces  is  continuous,  as  is  the  evaluation  function. 

In  the  classical  differential  calculus,  differentials  of  a  small  class  of  basic  functions  are  obtained 
ab  initio.  Differentials  of  functions  generated  from  the  basic  functions  by  generalized  composition 
are  obtained  via  the  chain  rule.  Differentiation  factors  out  of  composition,  and  the  identity  func¬ 
tion  is  its  own  differential  everywhere.  It  is  often  underappreciated  by  students  of  calculus  that 
familiar  rules  follow  from  the  chain  rule  and  details  of  the  differentials  of  various  basic  functions. 
For  example,  the  product  rule  follows  from  the  chain  rule  and  the  differential  of  the  multiplica¬ 
tion  operation.  Also,  proofs  of  the  chain  rule  in  elementary  calculus  texts  often  adopt  unneces¬ 
sary  hypotheses  of  continuous  differentiability.  These  proofs  do,  correctly,  use  the  continuity  of 
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multiplication,  but  this  actually  is  a  projection  of  the  continuity  of  the  composition  operation  on 
differentials  in  CONV. 

Category-theoretic  characterizations  of  heterotic  models  of  computation  were  introduced  by 
Stepney  et.al.  in  [40].  Heterotic  models  of  computation  combine  computational  models  such 
as  classical/quantum,  digital/analog,  synchronous/asynchronous,  imperative/functional/relational, 
etc.  to  obtain  increased  computational  power,  both  practically  and  theoretically.  In  the  project  we 
had  been  developing  heterotic  dynamical  systems  consistent  with  Stepney  et.al.’s  notions,  but  had 
been  unaware  of  their  work  and  terminology  until  after  their  work  appeared.  These  systems  involve 
the  evolution  of  closely  coupled  variables  of  broadly  varying  type,  some  discrete-valued,  some 
continuous-valued,  and  some  with  values  that  are  not  quantitative.  The  definition,  representation 
and  manipulation  of  differentials  is  non-obvious  in  such  circumstances. 

There  is  another  important  observation  that  is  obscured  in  the  elementary  calculus  of  real¬ 
valued  functions  on  the  real  numbers  (elementary  calculus  in  one  dimension):  derivatives  of  func¬ 
tions  that  are  differentiable  on,  say,  an  interval,  are  functions  that  are  differential  valued.  (This  is 
already  observable  in  multi-dimensional  elementary  calculus.)  The  function  space  of  differentials 
arising  in  this  context  is  the  space  of  linear  functions  from  the  reals  to  the  reals  as  a  subspace  of 
the  convergence  space  of  continuous  functions  from  the  reals  to  the  reals  in  CONV.  This  space  of 
differentials  is  homeomorphic  to  the  space  of  real  numbers  equipped  with  the  standard.  Euclidean , 
topology.  Using  the  homeomorphism,  one  then  identifies  the  differentials  with  real  numbers  and 
thereby  obtains  real-valued  derivatives.  Much  of  this  is  invisible  in  standard  calculus  texts,  even  at 
fairly  advanced  levels.  It  is  important  to  note  that  in  more  general  contexts  the  homeomorphisms 
through  which  such  identifications  are  made  will  not  always  exist.  In  particular,  notions  of  contin¬ 
uous  differentiability  then  require  that  derivatives  be  differential- valued  and  require  that  something 
like  topological  structure  (in  our  case,  convergence  structure)  on  the  function  spaces  that  arise  be 
taken  into  consideration. 

If  it  is  insisted  upon  that  differentials  be  linear  functions,  then  several  options  present  them¬ 
selves.  Most  stringently,  continuous  ring-module  structures  can  be  extracted  from  the  symmetries 
in  the  convergence  structures  on  the  spaces  involved.  Less  stringently,  functions  in  the  role  of 
linear  functions  are  determined  by  the  continuous  homormorphisms  among  the  translation  groups 
(< cf  section  on  the  spaces  involved.  When  the  translation  groups  are  both  full  and  Abelian  we 
obtain  continuous  ring  modules.  However,  blind  insistence  upon  linearity  is  misguided.  Even  in 
Euclidean  spaces,  one  obtains  an  interesting  differential  calculus  by  taking  the  continuous  homo¬ 
geneous  functions  of  order  1  (i.e.  f(ax)  =  af(x ))  as  the  differentials.  One  even  obtains  a  more 
general  but  less  interesting  calculus  by  dropping  the  continuity  requirement  on  the  differentials, 
but  preserving  continuity  of  composition. 

Over  time,  a  number  of  researchers  have  sought  to  generalize  differentiability  to  spaces  where 
the  generalization  is  non-obvious.  Some  of  the  more  serious  and  sophisticated  results  in  this  di¬ 
rection  have  employed  one  or  another  restriction  of  the  notion  of  convergence  space,  often  near 
to  pre-topological  spaces,  or  else  stayed  within  TOP  [7,  8,  11,  12,  22,  23,  28,  31,  34,  35].  These 
explorations  assumed  the  existence  of  additional  structure  characterizing  linearity.  [31]  recognized 
the  importance  of  Cartesian-closure  for  obtaining  a  robust  chain-rule. 

3.2  Convergence  spaces:  an  overview 

A  filter  on  a  set  X  is  a  collection  of  subsets  of  X  closed  under  finite  intersection  and  reverse 
inclusion.  A  proper  filter  is  a  filter  which  does  not  contain  the  empty  set  as  a  member.  Let  $(3f) 


Approved  for  Public  Release;  Distribution  Unlimited. 

6 


denote  the  set  of  all  filters  on  X.  For  a  subset  A  of  X,  {B  \  A  C  B  C  A"}  is  a  member  of  <1>(A"). 
We  denote  this  filter  by  [A].  In  the  special  case  where  A  is  a  singleton  {x}  we  denote  [A]  by  [x] 
and  call  this  the  point  filter  at  x. 

Definition  3.1.  A  convergence  structure  on  X  is  a  relation  |  (read  as  “converges  to”)  between 
members  of  <3>(A")  and  members  of  X  such  that  for  each  x  G  X:  (1)  [x]  converges  to  x,  and  (2)  the 
set  of  filters  converging  to  x  is  closed  under  reverse  inclusion.  A  pair  (X, /)  consisting  of  a  set  X 
and  a  convergence  structure  |  on  X  is  called  a  convergence  space. 

A  function  /  :  X  — >  Y  where  A"  and  Y  are  sets,  induces  functions  /  :  2X  — >  2)  and 

/  :  $(X)  — y  <3>(Y).  /  is  defined  by  f(A)  =  {/(a)  |  a  G  A},  which  we  call  the  f -image  of  A 
w.r.t  /.  For  F  G  T>(X)  note  that  the  collection  of  all  supersets  of  /-images  of  members  of  F  forms 

a  filter  which  we  call  f(F).  Hereafter  we  overload  notation  and  drop  the  "and  "annotations. 

When  convenient,  we  will  refer  to  a  convergence  space  (X,  /  )  by  its  carrier,  A". 

Definition  3.2.  Let  /  :  X  — >■  Y  where  X  and  Y  are  convergence  spaces,  and  let  x()  G  A".  /  is 
continuous  at  x0  iff  for  each  F  G  $(X),  if  F  j.  x0  in  A",  then  f(F )  j,  f(x0 )  in  Y.  f  is  continuous 
iff  /  is  continuous  at  every  point  of  X. 

Continuity  can  be  characterized  in  terms  of  filter  members,  which  play  a  role  analogous  to  the 
role  played  by  neighborhoods,  as  supersets  of  open  sets,  in  topological  spaces. 

Proposition  3.1.  Let  /  :  X  — >  Y  where  X  and  Y  are  convergence  spaces,  and  let  x0  be  a  point  of 
X.  f  is  continuous  at  x0  iff  for  every  filter  F  converging  to  x0  in  X,  there  is  a  filter  Q  converging 
to  f(x o)  in  Y  such  that  (W  G  Q)  (3 U  G  F)  [f(U)  C  V],  | 

Definition  3.3.  A  homeomorphism  between  two  convergence  spaces  is  a  continuous  bijection 
whose  inverse  is  continuous. 

Definition  3.4.  Let  a;  be  a  point  of  a  convergence  space  A",  and  let  U  be  a  subset  of  X.  U  is  said  to 
be  a  neighborhood  of  x  iff  U  belongs  to  every  filter  converging  to  x. 

Definition  3.5.  [19]  A  convergence  space  (A",  .|)  is  said  to  be  a  pretopological  space  if  and  only  if 
■I  is  a  pretopology ,  i.e.  for  each  x  G  A",  the  collection  of  all  neighborhoods  of  x  converges  to  x. 

Proposition  3.2.  Let  /  :  X  — >■  Y  where  A"  and  Y  are  pretopological  spaces,  and  let  x0  G  A".  /  is 
continuous  at  x0  iff  for  every  neighborhood  V  of  f(x0),  there  is  a  neighborhood  U  of  x0  such  that 
f(U)  CL  I 

It  is  evident  that  every  topological  space  is  a  pretopological  space  (cf.  [19,  18,  29]),  and  that 
the  convergence  space  notion  of  continuity  and  the  topological  space  notion  of  continuity  coincide 
for  topological  spaces.  As  indicated  in  the  introduction,  the  spaces  and  functions  of  interest  to  us 
are  naturally  organized  as  categories.  The  main  categories  of  interest  in  this  paper  are: 

CONV  the  category  of  convergence  spaces  and  continuous  functions 
PreTOP  the  category  of  pretopological  spaces  and  continuous  functions 
TOP  the  category  of  topological  spaces  and  continuous  functions 

RDIGRAPH  the  category  of  reflexive  digraphs  (i.e.  directed  graphs 

with  a  loop  at  each  vertex)  and  edge-preserving  functions 
PostDISCRETE 

the  full  subcategory  of  CONV  whose  objects  are  the 
post  discrete  convergence  spaces 
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TOP  is  a  full  subcategory  of  PreTOP,  which,  in  turn,  is  a  full  subcategory  of  CONV.  Both  of 
these  full  inclusions  are  reflective,  via  induced  pretopology  and  induced  topology  operations,  re¬ 
spectively.  RDIGRAPH  is  isomorphic  to  PostDISCRETE,  which,  in  turn,  embeds  into  PreTOP. 1 

The  reflection  functor  PreT  from  CONV  to  PreTOP  can  be  realized  by  letting  the  carrier  of 
PreT(X)  be  the  carrier  of  X,  and  and  letting  a  filter  T  converge  to  a  point  x  in  PreT(X)  iff  the 
collection  of  all  neighborhoods  of  x  in  X  is  a  subcollection  of  T . 

Similarly,  the  reflection  functor  T  from  PreTOP  to  TOP  can  be  realized  by  letting  the  carrier 
of  T(X)  be  the  carrier  of  X,  and  defining  the  topology  on  T(X)  as  {U  C  X  \U  is  a  neighborhood 
in  X  of  each  point  of  U} 

3.3  Reflexive  digraphs  as  convergence  spaces 

RDIGRAPH  can  be  embedded,  in  more  than  one  way,  as  a  full  subcategory  of  CONV. 

Definition  3.6.  A  convergence  space  X  will  be  said  to  be  postdiscrete  if  and  only  if  every  conver¬ 
gent  proper  filter  is  a  point  filter. 

Proposition  3.3.  The  postdiscrete  pretopological  spaces  are  precisely  the  discrete  topological 
spaces.  | 

Definition  3.7.  Let  (V,  E)  be  a  reflexive  digraph.  Induce  a  convergence  structure  on  V  by  letting  a 
proper  filter  T  converge  to  a  vertex  x  iff  J7  =  [y\  for  some  vertex  y  with  an  edge  in  E  from  x  to  y. 

It  is  readily  verified  that  if  (Vi,  E\ )  and  (V2,  E2)  are  reflexive  digraphs,  then  a  function  /  :  Vi  — y 
V2  is  continuous  (with  respect  to  the  induced  convergence  structures  on  Vi  and  V2)  iff,  for  all  edges 
( x ,  y)  in  Ei,  the  edge  (f(x),  f(y))  is  present  in  E2. 

Proposition  3.4.  The  construction  in  Definition  3.7  embeds  RDIGRAPH  as  a  full  subcategory  of 
CONV,  namely  the  full  subcategory  whose  objects  are  the  postdiscrete  spaces. 

This  embedding  is  coreflective.  The  coreflection  functor  RDIGR  APH  :  CONV  — »  RDIGRAPH 
can  be  obtained  letting  the  vertices  of  RDIGRAPH (X)  be  the  members  of  X,  and  letting  an  or¬ 
dered  pair  (x,  y)  be  an  edge  of  RDIGRAPH(X)  iff  [y]  j,  x  in  X.  | 

Alternatively,  RDIGRAPH  can  be  embedded  as  a  full  subcategory  of  PreTOP,  and  thence  as  a 
full  subcategory  of  CONV  [41,  39],  by  letting  a  filter  T  converge  to  a  vertex  x  iff  {y  |  (x,  y)  G  E} 
is  a  member  of  T . 

In  general,  this  embedding  of  RDIGRAPH  into  PreTOP  imposes  a  weaker  convergence  struc¬ 
ture  on  reflexive  digraphs  than  the  embedding  in  Definition  3.7. 

Proposition  3.5.  The  embedding  of  RDIGRAPH  into  PreTOP  [41,  39]  is  the  composite  of  the 
embedding  in  Definition  3.7  of  RDIGRAPH  into  CONV  with  the  reflection  functor  from  CONV 
to  PreTOP,  and  embeds  RDIGRAPH  as  a  full,  coreflective  subcategory  of  PreTOP,  and  thence 
as  a  full,  coreflective  subcategory  of  CONV. 

The  coreflection  functor  from  PreTOP  to  RDIGRAPH  (and  the  coreflection  functor  from 
CONV  to  RDIGRAPH  via  PreTOP)  can  be  obtained  in  precisely  the  same  way  as  in  Proposi¬ 
tion  3.4.  | 

'The  embedding  is  the  restriction  of  the  induced  pretopology  reflection  from  CONV  to  PostDISCRETE 
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The  reflexive  digraphs  whose  induced  pretopologies  are  topological  are  precisely  those  in 
which  the  underlying  binary  relation  is  transitive  as  well  as  reflexive.  [41,  39] 

3.4  Function  spaces 

Unlike  TOP  and  PreTOP,  CONV  is  a  Cartesian  closed  category  ([33,  27]): 

Definition  3.8.  [27]  Let  X  and  Y  be  convergence  spaces. 

The  function  space  Yx  is  the  set  of  all  continuous  functions  from  X  to  Y,  equipped  with  the 
convergence  structure  j.  defined  as  follows: 

For  each  'H  G  <f>(FA  )  and  each  f0  G  Yx,  let  H  |  ./0  if  and  only  if 
for  each  x0  G  X  and  each  T  j.  x0,  {  {  f(x)  \f  G  F[, x  E  F}\F[  e  PI,  FeF}  is  a  base  for  a 
filter  which  converges  to  f(x0)  in  Y. 

Definition  3.8  is  essentially  the  same  as  the  definition  of  function  graphs  in  [26].  This  should 
not  be  very  surprising:  RDIGRAPH  is  a  Cartesian  closed  category  (cf.  [3]). 

Proposition  3.6.  [11,  13]  Let  A"  and  Y  be  convergence  spaces.  Then  the  evaluation  function  from 
Yx  x  X  to  Y  is  continuous.  | 

The  following  is  a  very  useful  consequence  of  Proposition  3.6. 

Corollary  3.1.  Let  X,  Y ,  and  Z  be  convergence  spaces.  Then  the  composition  function  from 
Z't  x  yx  to  Zx  is  continuous.  | 

3.5  (Appendix  to  section  3:  Reflective  and  coreflective  subcategories) 

The  content  of  this  section,  including  the  second  example,  is  largely  taken  from  [33].  The  first 
example  is  taken  from  [25].  Also  see  [6]  and  [3]. 

Throughout  this  section  let  A  be  a  subcategory  of  a  category  B. 

A  is  said  to  be  a  reflective  subcategory  of  B  iff  there  is  a  reflection  functor  from  B  to  A,  i.e.  a 
functor  F  :  B  — *  A  such  that  each  object  B  of  B  is  equipped  with  a  reflector,  i.e.  an  arrow 
tb  :  B  — >■  FB  in  B  such  that,  for  each  object  A  of  A  and  each  arrow  /  :  B  — >  A  in  B,  there  is 
a  unique  arrow  /  :  FB  — )■  A  in  A  such  that 


B  —A  pB 


C 


commutes. 

If  A  is  a  full,  reflective  subcategory  of  B  with  reflection  functor  F.  then,  for  each  object  A  of  A, 
A  and  FA  are  isomorphic. 

For  instance,  the  category  of  compact  Hausdorff  spaces  and  continuous  functions  is  a  full,  reflective 
subcategory  of  Top.  The  reflection  functor  fl  is  the  well-known  Stone-Cech  functor.  The  reflector 
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from  an  arbitrary  topological  space  X  to  f3X  is  an  embedding  if  and  only  if  X  is  completely  regular 
and  Hausdorff. 


Dually,  A  is  said  to  be  a  coreflective  subcategory  of  B  iff  there  is  a  coreflection  functor  from  B  to 
A,  i.e.  a  functor  F  :  B  — >  A  such  that  each  object  B  of  B  is  equipped  with  a  coreflector,  i.e.  an 
arrow  lb  '■  FB  — *  B  in  B  such  that,  for  each  object  A  of  A  and  each  arrow  /  :  A  — *  B  in  B, 
there  is  a  unique  arrow  /  :  A  — »  FB  in  A  such  that 

FB,  >  B 
f 
C 


commutes. 

If  ^4  is  a  full,  coreflective  subcategory  of  B  with  coreflection  functor  F.  then,  for  each  object  A  of 
A,  A  and  FA  are  isomorphic. 

For  instance,  the  category  of  Abelian  torsion  groups  and  homomorphisms  is  a  full,  coreflective 
subcategory  of  the  category  of  Abelian  groups  and  homomorphisms. 

4  RESULTS  AND  DISCUSSION 

4.1  Translation  groups  and  homogeneous  convergence  spaces 

Definition  4. 1.  An  automorphism  of  a  convergence  space  X  is  a  homeomorphism  /  :  X  — *  X. 

Definition  4.2.  A  translation  group  on  a  convergence  space  X  is  a  group  T  of  automorphisms  of 
X  such  that,  for  each  pair  of  points  p  and  q  of  X,  there  is  at  most  one  member  of  T  which  maps  p 
to  q.  In  general,  we  will  denote  this  unique  member  of  T  (if  it  exists)  by  (q  —  p). 

Notation:  The  group  operation  of  a  translation  group  T  on  a  convergence  space  X  will  be  written 
additively,  whether  or  not  T  is  Abelian.  Furthermore,  for  all  r  €  T  and  all  x  €  A",  we  will  write 

t{x)  as  x  +  r 

In  this  notation,  the  requirement  that  the  translation  (q  —  p)  (if  it  exists)  maps  p  to  q  becomes 
the  familiar  requirement  that  if  (q  —  p)  exists,  then  p  f  (q  —  p)  —  q 

A  full  translation  group  on  a  convergence  space  X  is  a  translation  group  on  X  which  contains 
a  translation  (q  —  p)  for  each  pair  of  points  p  and  q. 

Proposition  4.1. 

i.  Every  convergence  space  X  can  be  embedded  as  a  subspace  of  a  convergence  space  HX 
which  has  a  full  translation  group. 

ii.  X  and  HX  have  the  same  cardinality  if  and  only  if  the  cardinality  of  X  is  either  zero  or 
infinite. 

iii.  The  embedding  of  X  into  HX  is  onto  HX  if  and  only  if  X  is  empty. 
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iv.  If  X  and  Y  are  arbitrary  convergence  spaces  then  every  continuous  function  /  :  X  — »  Y 
can  be  be  extended  to  a  continuous  function  Hf  :  HX  — >  HY . 

v.  If  /  is  a  homeomorphism,  then  so  is  Hf. 


An  immediate  consequence  of  (ii)  in  Proposition  4. 1  is  that,  if  X  is  a  finite  space  with  (or 
without)  a  full  translation  group,  then  HX  cannot  be  homomeomorphic  to  X.  It  should  also  be 
noted  that,  given  a  particular  /,  the  continuous  extention  Hf  need  not  be  unique. 

Definition  4.3.  A  convergence  space  X  is  homogeneous  iff  for  each  pair  of  points  x\  and  x2  of  X, 
there  is  an  automorphism  of  X  which  maps  X\  to  x-2 

Obserx’ation  4.1.  A  convergence  space  which  has  a  full  translation  group  must  be  homogeneous. 
Furthermore,  a  full  translation  group  on  a  nonempty  convergence  space  X  must  have  the  same 
cardinality  as  X.  | 

4.2  Differential  Calculi 

Definition  1.  A  differential  calculus  is  a  category  V  in  which 

i.  every  object  of  V  is  a  triple  X  =  (X,  0,  T)  such  that  X  is  a  convergence  space,  0  is  a  point 
of  X  (called  the  origin  of  X),  and  T  is  a  full  translation  group  on  X. 

ii.  every  arrow  in  V  from  an  object  (. X ,  Ox,  Tx)  to  an  object  (Y.  0y ,  TY )  is  a  continuous  function 
from  X  to  Y  which  maps  Ox  to  0y 

iii.  composition  of  arrows  in  V  is  function  composition. 

iv.  for  every  object  X  =  (. X ,  Ox,  Tx),  the  identity  function  on  X  is  an  arrow  in  V  from  X  to  X 

v.  for  each  pair  of  objects  X  =  (. X ,  Ox,  Tx )  and  y  =  (Y,  0y ,  TY),  the  constant  function  map¬ 
ping  every  point  of  X  to  0y  is  an  arrow  in  V  from  X  to  y 


In  view  of  Proposition  4.1,  the  requirement  that  each  object  have  a  full  translation  group  is  not 
unduly  restrictive. 

4.3  Differentials 

Let  a  G  A  C  X  and  let  B  C  Y,  where  X  =  (A",  Ox,  Tx)  and  y  =  (Y,  0y,  TY)  are  objects  of  a 
differential  calculus  V.  Let  /  :  A  — »  B  be  an  arbitrary  function. 

Let  L  e  T>(X,  y),  where  T>(X,  y)  is  the  set  of  all  arrows  in  V  from  X  to  y,  equipped  with  the 
subspace  convergence  structure  inherited  from  the  function  space  Y'x  in  CONV. 

Definition  2.  L  is  a  differential  of  /  at  a  iff 

for  every  T  f  a  in  A,  there  is  some  H  f  L  in  V(X,  y)  such  that 

i.  H  C  [L],  and 
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ii.  for  every  H  e  EL,  there  is  some  F  e  F  such  that 

for  every  point  x  £  F,  there  is  at  least  one  function  t  e  H  such  that 

t(x  +  (0X  -  a))  =  f(x)  +  (Oy  -  /(a)) 

(In  Definition  2,  (/(a)  —  Oy)  o  t  o  (Ox  —  o)  is  an  extrapolant  of  /  through  (a,  /(a))  and  (xj(x)).) 

Definition  3.  A  function  from  /I  to  B  is  differentiable  (respectively,  uniquely  differentiable )  at  a 
point  a  iff  it  has  at  least  one  (respectively,  precisely  one)  differential  at  a. 

A  function  from  /I  to  B  is  differentiable  (respectively,  uniquely  differentiable)  iff  it  is  differ¬ 
entiable  (respectively,  uniquely  differentiable)  at  each  point  of  A. 

4.4  The  Chain  Rule 

As  we  indicated  in  the  introduction,  the  chain  rule  is  among  the  central  roles  in  our  differential 
calculi.  In  elementary  calculus  for  example,  the  product  rule  follows  from  the  chain  rule  after 
obtaining  the  differential  of  the  multiplication  operation. 

Example  4.1.  Expressed  in  terms  of  differentials,  the  product  rule  for  real- valued  functions  of  a 
real  variable  reduces  to  matrix  multiplication  (i.e.  composition  of  linear  functions). 

Dx (mult  O  (/,  g))  =  DUtg)(x) mult  o  Dx(f,  g ) 

=  A/(*),s(*))mu It  °  (Dxf,  Dxd) 

=  \g(x)f(x))\Dff 
=  g{x)Dxf  +  f(x)  Dxg 

Returning  to  our  more  general  setting,  let  a  e  A  C  X,  let  B  C  Y ,  and  let  C  C  Z,  where 
X  =  (X.  Ox,  Tx ) ,  y  =  (Y.  Oy,  Ty),  and  Z  =  ( Z ,  0Z,  Tz)  are  objects  of  a  differential  calculus  V. 
Let  /  :  A  — >  B  and  g  :  B  — >  C  be  arbitrary  functions. 

Let  K  :  X  — >  y  and  L  :  y  — >  Z  be  arrows  of  V 

Theorem  4. 1.  (Chain  Rule)  Suppose  that  /  is  continuous  at  a.  Also  suppose  that  K  is  a  differential 
of  /  at  a,  and  L  is  a  differential  of  g  at  f(a). 

Then  L  o  K  is  a  differential  of  g  o  /'  at  a. 

proof:  Let  F  be  a  filter  converging  to  a  in  X.  Since  K  is  a  differential  of  /  at  a,  there  is  some 
Q  l  K  in  V(X,  y)  such  that  Q  C  [K]  and,  for  every  G  e  Q.  there  is  some  F{  G  e  F  such  that  for 
each  point  x  G  Fi  g  there  is  some  function  sq, x  G  G  such  that 

SG,x(x  +  (°X  -  a))  =  f(x)  +  (Oy  -  /(a))  (1) 

On  the  other  hand,  since  /  is  continuous  at  a,  we  have  /(F)  j,  /(a)  in  B.  Since  L  is  a 
differential  of  g  at  /(a),  there  is  some  filter  EL  l  L  in  V(y.  Z)  such  that  EL  C  [L]  and,  for  every 
H  e  EL,  there  is  some  Nh  €  /(F)  such  that  for  each  point  y  e  jV#,  there  is  some  function 
tn.y  €  V  such  that 

*/f,i/(j/  +  (°v  -  /(a)))  =  ^(/(®))  +  (0 z  ~  ff(f(a))) 
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(2) 


Consider  such  a  set  NH.  By  definition,  NH  E  f{F),  i.e.  there  is  some  F2H  G  F  such  that 

f{F2,H)  C  iVH 

By  (2),  for  each  point  x  G  we  have 

+  (Oy  -  /(a)))  =  c/(/(x))  +  (0 z  -  g(f(a )))  (3) 

Next,  note  that  {  {  ft,2  °  /ii  |  hi  E  G,h2  G  H  }  |  G  G  (/,  H  G  "H  }  is  a  basis  for  a  filter  J-  on 
£>(*,  Z),  and  that  J  C  [L  o  K\. 

By  joint  continuity  of  composition  (Corollary  3.1),  J  j.  L  o  L  in  V(X,  Z). 

Let  J  be  an  arbitrary  member  of  J .  There  exist  G  G  Q  and  II  G  'H  such  that 

{  h2  o  h\  |  hi  G  G,  h2  G  H  }  C  J 

Let  F  =  F\  g  fl  F2ih-  Then  F  G  T .  For  each  point  x  G  F,  we  have  G  G  and  G  H, 

and  therefore  tHj(x )  °  G  J.  Furthermore,  by  (1)  and  (3), 

tH,f(x)(sG,x(x  +  (°x  -  «)))  =  tHj(x)(f(x)  +  (Oy  -  /(a)))  =  ^(/(x))  +  (0Z  -  g(f(a )))  (4) 

Thus,  (g{f(a))  —  0z)°tH,f(x)°SG,xo{0x—a)  is  the  required  extrapolant  of  gof  through  ( a,g(f(a ))) 
and  (x,g(f(x))). 

Since  the  selection  of  x  is  arbitrary  (once  G  and  II  have  been  chosen),  L  o  K  is  indeed  a 
differential  of  g  o  /  at  a.  | 

4.5  Examples  of  Differential  Calculi 

In  all  examples,  let  R  be  the  real  line,  and  let  N  be  the  set  of  all  natural  numbers. 

Example  4.2.  The  classical  Newton-Leibniz  differential  calculus  of  real  variables 

The  objects  of  this  differential  calculus  are  the  spaces  Rn  (n  G  N),  equipped  with  their  usual 
(Euclidean)  convergence  structures,  their  usual  origins,  and  their  usual  translation  groups. 

The  arrows  of  this  calculus  are  the  R-linear  functions. 

In  this  calculus,  differentiability  and  unique  differentiability  are  equivalent,  and  a  function  / 
has  a  differential  at  a  point  p  iff  /  is  differentiable  (in  the  usual  sense)  at  p. 

Example  4.3.  A  directional  calculus  of  real  variables 

This  is  an  extension  of  the  Newton-Leibniz  calculus  of  real  variables.  Again,  the  carriers 
of  objects  are  the  sets  R"  (n  G  N),  equipped  with  their  usual  origins  and  their  usual  translation 
groups. 

However,  each  R"  will  be  equipped  with  two  convergence  structures,  and  thus  will  be  be  the 
carrier  of  two  distinct  objects. 

One  of  these  convergence  structures  is  the  familiar  Euclidean  structure. 

The  other  convergence  structure  which  we  impose  on  Rn  is  defined  as  follows: 

A  filter  F  will  be  said  to  converge  to  a  point  p  iff  there  is  some  unit  vector  q  such  that 

{p  +  aq\a  E  H,  \a\  <  e}  E  F 

for  every  real  number  e  >  0.  We  call  this  convergence  structure  the  directional  con¬ 
vergence  structure  on  R". 
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For  n  —  1,  the  directional  convergence  structure  coincides  with  the  Euclidean  structure.  For 
n  >  1,  this  directional  convergence  structure  is  stronger  than  the  Euclidean  structure,  and  is  not 
pretopological. 

The  arrows  of  this  calculus  are  the  continuous  R-homogeneous  functions  of  degree  one. 

Again,  differentiabilty  and  unique  differentiabilty  are  equivalent.  Observe  that  a  function  / 
from  Rm  with  its  directional  convergence  structure  to  Rn  with  its  Euclidean  convergence  structure 
has  a  differential  at  a  point  p  iff  /  has  directional  derivatives  in  all  directions  at  p. 


Example  4.4.  The  one-sided  calculus  of  a  single  function  of  a  real  variable 

The  only  object  of  this  differential  calculus  is  the  real  line  R,  equipped  with  its  usual  (Eu¬ 
clidean)  convergence  structure,  its  usual  origin,  and  its  usual  translation  group. 

Define  the  arrows  from  R  to  itself  to  be  the  set  of  all  functions  L  :  R  — >■  R  such  that 


for  all  x,  y  G  R  and  every  a  >  0 

L(a(x  +  y))  =  aL(x)  +  aL(y) 


In  this  calculus,  a  function  has  a  differential  at  a  point  p  iff  it  has  both  a  left-sided  derivative 
at  p  and  a  right-sided  derivative  at  p.  As  in  the  previous  examples,  differentiability  and  unique 
differentiability  are  equivalent. 

Example  4.5.  The  Reed-Akers-Lee-Vichniac  calculus  of  Boolean  partial  derivatives 

([38,  4,32].) 

Let  B  be  the  two-element  field,  i.e.  the  field  of  integers  modulo  two. 

Let  m  be  a  natural  number,  and  let  /  be  an  arbitrary  function  from  Bm  to  B.  Following  [38],  [4], 
and  [32],  define  the  /'1h  partial  derivative  of  /  at  a  point  x  as  the  difference  quotient 


,  f  (xi,  X2-,  ■  ■  ■  i  Xm )  f  (3^1,  X2 ,  ‘  ■  1  1,  ~ 'Xj,  •  •  •  ;  Xm ) 

Tj  \x 

OXi  Xi  —  -IX  i 


Observe  that  this  is  always  well-defined.  Also  observe  that,  since  the  arithmetic  is  modulo  two, 
this  definition  simplifies  to 


df_ 

d  Xi 


X 


f  (*^1 ,  X2 1  •  •  •  1  Xm )  T  f  (x  1 ,  X2 1  •  ■  ■  1  Xi— 1 ,  '3?i ,  Xi- |_i ,  Xi-\- 2 1  •  •  •  1  %mri) 


In  other  words,  the  ith  partial  derivative  of  /  at  (xi,x2,  ■  ■  ■  xm)  is  1  if  changing  the  value  of  xt 
changes  the  value  of  f(x i,x2,  ■  ■  ■ ,  xm),  and  is  0  otherwise. 

Following  [32]  and  [43],  if  m  and  n  are  natural  numbers,  define  the  Jacobean  matrix  of  an  arbitrary 
function  F  from  Bm  to  B"  as  the  matrix  of  partial  derivatives  of  the  component  functions  of  F. 
Note  that  this  too  is  always  well-defined.  2 

Higher-order  partial  derivatives,  mixed  partials,  etc.,  are  defined  in  [9]  in  the  obvious  way. 

Many  of  the  rules  of  the  classical  differential  calculus  of  real  variables,  such  as  the  product  rule, 
also  hold  for  Boolean  derivatives  [9].  However,  it  is  notorious  that  the  conclusion  of  the  chain  rule 
does  not  always  hold,  if  partial  derivatives  are  involved  [43].  We  can  explain  this. 

2In  [32],  the  Jacobean  matrix  of  F  is  called  the  differential  of  F. 
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For  each  natural  number  n,  take  the  members  of  Bn  as  the  vertices  of  a  digraph,  with  an  edge  from 
a  vector  a  =  (ai,  a2, . . . ,  an)  to  a  vector  b  =  (b\,  b2, . . . ,  bn)  iff  the  Hamming  distance  between  a 
and  b  is  at  most  1,  i.e.  there  is  at  most  one  i  such  that  a(i )  ^  b(i). 

In  other  words,  take  the  members  of  Bn  as  the  vertices  of  a  hypercube  with  two-way  edges,  and 
with  a  loop  at  each  vertex.  Equip  Bn  with  the  postdiscrete  convergence  structure  induced  by  this 
digraph  (cf.  Definition  3.7). 

In  short,  the  proper  filters  which  converge  to  a  bit  vector  a  are  precisely  the  point  filters  [6]  such 
that  the  bit  vector  b  differs  from  a  in  at  most  one  bit. 

Take  these  spaces  B"  as  the  carriers  of  the  objects  of  a  differential  calculus.  For  each  n,  take  the 
zero  vector  of  length  n  as  the  origin  of  Bn. 

The  group  generated  by  the  flips  hi,h2,...,  hn  is  taken  as  the  translation  group  of  Bn,  where  (as 
one  would  expect)  h^(b)  is  obtained  from  bit  vector  b  by  changing  the  kth  bit  of  b  (and  leaving 
every  other  bit  unchanged). 

For  each  m  and  n,  take  the  continuous  B-linear  functions  from  Bm  to  B"  as  the  differentials  from 
Bm  to  B"  . 

In  particular,  there  are  precisely  two  differentials  between  B  and  itself,  namely,  the  constant  func¬ 
tion  which  returns  zero,  and  the  identity  function. 

In  this  differential  calculus,  as  in  the  previous  two  examples,  differentiability  and  unique  differen¬ 
tiability  are  equivalent. 

It  is  readily  verified  that  if  a  function  F  from  Bm  to  Bn  has  a  differential  L  at  a  bit  vector  x,  then 
L  is  the  B-linear  function  whose  matrix  (with  respect  to  the  standard  bases  of  Bm  and  Bn)  is  the 
Jacobean  matrix  of  F. 


Function 

Differential  at  0 

Differential  at  1 

constant  with  value  0 

constant  with  value  0 

constant  with  value  0 

constant  with  value  1 

constant  with  value  0 

constant  with  value  0 

identity 

identity 

identity 

— 1 

identity 

identity 

Figure  1:  Unary  Boolean  operators  and  their  Reed-Akers-Fee-Vichniac  hypercube  differentials. 
The  Reed-Akers-Fee-Vichniac  linear  differentials  of  unary  Boolean  operators  and  the  homoge¬ 
neous  differentials  of  unary  Boolean  operators  are  the  same  as  their  Reed-Akers-Fee-Vichniac 
hypercube  differentials. 


Every  function  between  B  and  itself  is  differentiable  everywhere. 
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Figure  2:  Boolean  connectives  and  their  hypercube  Reed- Akers-Lee- Vichniac  differentials. 


Likewise,  every  function  from  B2  to  B  is  everywhere  differentiable.  However,  not  every  function 
from  B  to  BJ  is  differentiable. 


To  see  that  Boolean  derivatives  do  not  always  obey  the  conclusion  of  the  chain  rule,  let  t  be  an 
independent  variable,  let  x  =  t,  let  y  =  t,  and  let  z  =  x\!  y.  Then 


since 


whereas 


dz 

dt 


dz  dx  dz  dy 
dx  dt  dy  dt 


dz 

dt 


dz  dx  dz  dy 
dx  dt  dy  dt 


0 


What  went  wrong? 

Answer:  Two  of  the  hypotheses  of  Theorem  4.1  do  not  hold  here.  Both  at  x  =  0  and  at  x  =  1, 


Dx(y  o  (id,  id))  =  Dxid  =  id 


but  we  cannot  assert  (indeed,  we  cannot  meaningfully  write) 

Dx(y  o  (id,  id))  =  (£>(*, *)V)  o  (Dx(id,id)) 

since  the  diagonal  embedding  (id,  id)  of  B  into  B2,  although  linear,  is  neither  continuous  nor 
differentiable  in  this  calculus. 
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We  can  remedy  this  failure  by  equipping  B2  with  the  product  convergence  structure  generated  by 
the  convergence  structure  on  B,  instead  of  the  hypercube  structure  used  in  Example  4.5. 

Example  4.6.  The  Reed-Akers-Lee-Vichniac  calculus  of  Boolean  total  derivatives 

([38,  4,32].) 

As  in  (Example  4.5),  take  the  finite-dimensional  vector  spaces  B"  as  the  carriers  of  objects,  take  the 
zero  vector  of  length  n  as  the  origin  of  B",  and  take  the  group  generated  by  the  flips  hi,  h2,  ■  ■  ■ ,  hn 
as  the  translation  group  of  Bn, 

As  before,  equip  the  two-element  field  B  with  the  postdiscrete  convergence  structure  induced  by 
taking  0  and  1  as  the  vertices  of  a  complete  digraph. 

However,  instead  of  equipping  Bn  with  a  hypercube  structure,  equip  Bn  with  the  product  con¬ 
vergence  structure  (which  coincides  with  the  postdiscrete  convergence  structure  induced  by  the 
complete  digraph  on  Bn). 

As  before,  for  each  rn  and  n,  take  the  continuous  B-linear  functions  from  Bm  to  Bn  as  the  differ¬ 
entials  from  Bm  to  B  \ 

In  general,  the  set  of  differentials  from  B"  to  B"  will  be  larger  than  before.  Every  function  between 
complete  digraphs  is  edge-preserving,  and  therefore  every  function  from  Bm  to  B"  is  continuous. 
In  short,  the  differentials  from  Bm  to  Bn  are  the  B-linear  functions. 

Since  the  convergence  structure  on  B  and  the  set  of  differentials  from  B  to  itself  are  the  same  as  in 
Example  4.5,  the  differentials  of  functions  from  B  to  itself  are  exactly  the  same  as  in  Example  4.5. 

Since  the  underlying  sets  of  objects  are  the  same  as  before,  and  there  are  more  differentials, 
one  might  perhaps  expect  differentiability  to  be  more  easily  satisfied  than  before,  but,  in  fact,  the 
reverse  is  true.  Differentiability  of  a  function  at  a  point  determined  by  the  function’s  local  behavior. 
Whereas,  before,  a  point  in  Bn  had  only  n  +  1  out-neighbors  (including  itself),  now  it  has  2”  out- 
neighbors.  In  short,  differentiability  of  a  funtion  at  a  point  now  depends  on  the  function’s  global 
behavior. 


We  leave  it  to  the  reader  to  verify  that  a  function  /  :  Bm  — >  Bn  is  differentiable  at  a  bit  vector  b 
if  and  only  if  /  is  an  affine  function,  i.e.  if  and  only  if  there  is  a  B-linear  function  L  :  Bm  — >  B" 
such  that 

f(x)=L(x)  +  f(  0)  (x  G  Bm) 

and  that,  if  such  an  L  exists,  then  L  is  the  unique  differential  of  /  everywhere  on  Bm. 

Again,  let  t  be  an  independent  variable,  let  x  —  t,  let  y  —  t,  and  let  z  =  x\l  y.  Since  partial 
derivatives  are  the  same  as  before,  it  is  still  true  that 


dz 

dt 


dz  dx  dz  dy 
dx  dt  dy  dt 


What  went  wrong  this  time? 

Answer:  Again,  there  is  a  hypothesis  of  Theorem  4.1  which  does  not  hold.  Both  at  x  =  0  and  at 

x  —  1, 

Dx(y  o  (id,  id))  =  Dxid  =  id 
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Figure  3:  Boolean  connectives  and  their  linear  Reed- Akers-Lee- Vichniac  differentials. 


but  again,  we  cannot  assert  (indeed,  cannot  meaningfully  write) 

Dx(y  o  (id,  id))  =  (£>(*, *)V)  o  (Dx(id,id)) 

The  diagonal  embedding  (id,  id)  of  B  into  BJ  is  continuous  (and,  in  fact,  is  a  differential)  in  this 
calculus,  but  V  is  nowhere  differentiable,  although  its  partials  exist  everywhere.  This  should  not 
be  astonishing;  analogous  phenomena  exists  in  the  classical  multivariable  calculus  of  Newton  and 
Leibniz. 

Nevertheless,  it  is  possible  to  define  a  Boolean  differential  calculus  with  a  reliable  chain  rule, 
provided  that  we  do  not  insist  that  differentials  be  linear. 

Example  4.7.  The  homogeneous  Boolean  differential  calculus 

Let  the  objects  be  precisely  the  same  as  in  Example  4.6. 

For  each  m  and  n,  define  the  differentials  from  Bm  to  B"  to  be  all  origin-preserving  functions 
from  Bm  to  Bn.  (Equivalently,  define  the  differentials  to  be  the  B-homogeneous  functions  of 
degree  one.)  As  noted  before,  this  is  compatible  with  our  definition  of  a  differential  calculus,  since 
every  function  between  complete  digraphs  is  edge-preserving. 

Since  the  convergence  structure  on  B  and  the  set  of  differentials  from  B  to  itself  are  unchanged, 
the  differentials  of  functions  from  B  to  itself  are  exactly  the  same  as  before. 
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Figure  4:  Boolean  connectives  and  their  homogeneous  differentials. 


The  situation  is  quite  different  for  functions  of  more  than  one  variable.  We  leave  it  to  the  reader 
to  verify  that  every  function  from  Bm  to  B  is  uniquely  differentiable. 

Since  all  functions  are  continuous  and  differential,  the  chain  rule  (Theorem  4.1)  is  always  applica¬ 
ble  in  the  homogeneous  Boolean  differential  calculus.  For  instance, 

id  =  -Do  id 

=  -D0(Vo  (id,  id)) 

=  (-D(o, o)V)  o  (D>0(id,id)) 

=  Vo  (id,  id) 

=  id 

Likewise, 

id  =  Di  id 

=  D>i(V  o  (id,  id)) 

=  (D(i;i)V)  o  (D>i(id,id)) 

=  A  o  (id,  id) 

=  id 
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4.6  Differential  Calculi  With  Nonhomogeneous  Objects 

In  the  preceding  development,  convergence  spaces  without  full  translation  groups  are  “second 
class  citizens”  in  the  sense  that  they  cannot  be  the  carriers  of  objects  of  a  differential  calculus.  A 
somewhat  more  general  (and  slightly  more  complicated)  concept  of  “differential  calculus”  permits 
all  convergence  spaces  to  be  carriers  of  objects. 

Obsen’ation  4.2.  Let  T  be  a  translation  group  on  a  convergence  space  X.  For  each  point  x  of  X, 
let  [x\T  be  the  T-  orbit  ofx,  i.e.  [x]T  =  {x  +  t  \  t  £  T  }. 

If  X  is  nonempty,  then  the  set  of  all  T  orbits  partitions  X  into  homogeneous  subspaces.  For 
each  T-orbit  [x]T,  the  restrictions  of  the  members  of  T  to  [x]T  form  a  full  translation  group  TLi  on 

NT- 

Each  T[x]  is  a  quotient  group  of  T.  | 

Definition  4.4.  A  system  of  origins  for  a  convergence  space  A"  with  respect  to  a  translation  group 
T  is  a  set  of  representatives  of  the  T-orbits,  i.e.,  a  subset  O  of  X  containing  precisely  one  member 
of  each  T -orbit. 

For  each  point  x  of  X,  let  ()x  be  the  unique  member  of  O  belonging  to  the  same  T-orbit  as  x. 

Definition  4.5.  Let  /  :  X  — *  Y  be  a  function  between  convergence  spaces.  Let  Ox  (Oy, 
respectively)  be  a  system  of  origins  for  X  with  respect  to  a  translation  group  S  (for  Y  with  respect 
to  a  translation  group  T,  respectively). 

i.  /  will  be  said  to  respect  orbits  iff,  for  each  pair  of  points  p  and  q  of  A",  if  p  and  q  lie  in  the 
same  S'-orbit,  then  fip)  and  f(q)  lie  in  the  same  T-orbit. 

ii.  /  will  be  said  to  be  preserve  origins  iff  f(Ox)  Y  Oy. 

Definition  4.6.  A  generalized  differential  calculus  is  a  category  'D  in  which 

i.  every  object  of  V  is  a  triple  X  =  (X,  T,  O )  such  that  A"  is  a  convergence  space,  T  is  a 
translation  group  on  X,  and  O  is  a  system  of  origins  for  X  with  respect  to  T. 

ii.  every  arrow  in  V  from  an  object  (X,  S,  Ox)  to  an  object  (Y.  T,  Oy)  is  a  continuous,  orbit- 
respecting,  origin-preserving  function  from  X  to  Y. 

iii.  composition  of  arrows  in  V  is  function  composition. 

iv.  for  every  object  X  =  (X,  T,  O),  the  identity  function  on  A"  is  an  arrow  in  V  from  X  to  X 

v.  for  each  pair  of  objects  X  =  (X,  S,  Ox)  and  y  =  (Y.  T,  Oy)  and  each  (  in  Oy ,  the  constant 
function  mapping  every  point  of  X  to  (  is  an  arrow  in  V  from  X  to  y 

A  differential  calculus  (in  the  sense  of  Definition  1)  is  essentially  the  same  notion  as  a  gen¬ 
eralized  differential  calculus  in  which  the  translation  group  of  every  object  is  a  full  translation 
group. 

At  the  opposite  extreme,  there  are  generalized  differential  calculii  in  which  the  translation 
group  of  every  object  is  trivial  (and  hence  all  orbits  are  singletons). 
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Example  4.8.  CONV  as  a  generalized  differential  calculus 

The  objects  of  the  trivial  generalized  differential  calculus  are  all  convergence  spaces,  equipped 
with  trivial  translation  groups.  The  arrows  from  an  object  X  to  an  object  Y  are  all  continuous 
functions  from  X  to  Y.  (Since  all  orbits  are  singletons,  every  function  is  orbit-respecting  and 
origin-preserving . ) 

Let  a  E  X  and  let  b  G  Y,  where  X  =  (X,  S,  Ox)  and  y  =  ( Y ,  T,  Oy)  are  objects  of  a 
generalized  differential  calculus  V.  Let  /  :  A  — »  B  be  an  arbitrary  function. 

Let  L  G  V(X,  y),  where,  again,  V(X,  y)  is  the  set  of  all  arrows  in  V  from  X  to  y,  equipped 
with  the  subspace  convergence  structure  inherited  from  the  function  space  Yx  in  CONV. 

Definition  4.7.  L  is  a  differential  of  /  at  a  iff 

for  every  T  |  a  in  X,  there  is  some  EL  |  L  in  'D(X .  >’)  such  that 

i.  TL  C  [{L}],  and 

ii.  for  every  H  G  EL,  there  is  some  F  G  T  such  that 

for  every  point  x  G  F,  there  is  at  least  one  function  t  G  H  such  that 

t(x  -  a)  =  f(x)  -  f(a) 


Note  that  y  —  a  =  ra  1(y). 

Differentiability  and  unique  differentiability  are  defined  exactly  as  before. 

Example  4.9.  The  classical  affine  differential  calculus  of  real  variables 

The  objects  of  this  generalized  differential  calculus  are  the  Euclidean  spaces,  equipped  with 
trivial  translation  groups.  The  arrows  from  Rm  to  R"  are  all  affine  functions  from  Rm  to  Rn 

As  in  the  classical  linear  differential  calculus,  a  function  /  has  a  differential  at  a  point  p  iff  / 
is  differentiable  (in  the  usual  sense)  at  p. 

Let  Epf  be  the  differential  of  /  at  p  in  the  classical  affine  differential  calculus  of  real  variables. 
I.e.,  Ep(f)  is  the  affine  function  which  best  approximates  /  in  arbitrarily  small  neighhborhoods  of 

p. 

Then  the  differential  of  /  at  p  in  the  classical  linear  differential  calculus  is  the  unique  linear 
function  which  can  be  obtained  from  Epf  by  composing  it  on  both  sides  with  translations  (in  the 
usual  sense),  i.e.  the  function  which  maps  each  point  q  to  f(p)  +  ( Epf)(q  —  p) 

4.7  The  Chain  Rule  for  Generalized  Differential  Calculii 

Let  a  G  X,  ,  where  X  =  (X,  R,  Ox),  V  =  (Y,S,Oy),  and  Z  =  (Z,T,Oz)  are  objects  of  a 
generalized  differential  calculus  V.  Let  /  :  X  — >  Y  and  g  :  Y  — »  Z  be  arbitrary  functions. 

Let  K  :  X  — »  y  and  L  :  y  — »  Z  be  arrows  of  V 

Theorem  4.2.  (Chain  Rule)  Suppose  that  /  is  continuous  at  a.  Also  suppose  that  K  is  a  differential 
of  /  at  a,  and  L  is  a  differential  of  g  at  f(a). 

Then  L  o  K  is  a  differential  of  g  o  /'  at  a. 

The  proof  of  Theorem  4.2  is  identical  to  the  proof  of  Theorem  4.1,  with  the  obvious  modifica¬ 
tions.  | 
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4.8  Heterotic  Dynamical  Systems:  Classical  and  Quantum 

Category-theoretic  characterizations  of  heterotic  models  of  computation  were  introduced  by  Step¬ 
ney  et  al.  in  [40].  As  remarked  above  in  section  3.1,  Heterotic  models  of  computation  combine 
computational  models  such  as  classical/quantum,  digital/analog,  synchronous/asynchronous,  im¬ 
perative/functional/relational,  etc.  to  obtain  increased  computational  power,  both  practically  and 
theoretically.  A  simple  but  highly  informative  classical/quantum  heterotic  model  of  computation 
is  represented  by  the  simple  sequential  imperative  quantum  programming  language  presented  by 
Abramsky  in  [2] . 

Abramsky’s  language,  in  its  simplest  form,  is  built  around  three  types:  nat,  bit  and  qubit.  The 
language’s  commands  are  of  six  kinds,  the  first  four  of  which  are:  assignment ,  noop ,  conditional 
and  loop  while.  The  ability  to  build  up  arithmetic  expressions  of  type  nat  and  boolean  expressions 
of  type  bit  is  presumed  and  intended  to  be  sufficient  for  having  a  “Turing  complete”  language. 
Commands  can  be  sequenced.  This  fragment  of  the  language  is  sequential  and  classical.  Finitely 
branching  probabilistic  nondeterminism  is  obtained  by  adding  the  command  colorfully  given  by 
cointoss  p  in  v,  where  p  is  an  explicit  real  number  in  the  closed  unit  interval  and  v  is  of  type  bit. 
Operationally,  when  executed,  the  cointoss  command  assigns  0  to  v  with  probability  p,  and  1  to 
v  with  probability  1  —  p.  A  stock  of  commands  consisting  of  1 -argument  and  2-argument  explicit 
unitary  qubit  operations  that  can  be  applied  to  variables  of  type  qubit  are  included.  Finally,  a 
quantum  measurement  command  is  included:  measure  q  in  v,  where  q  is  of  type  qubit  and  v  is 
of  type  bit. 

The  states  of  a  qubit  variable  are  wave  functions  of  the  classical  bit  values.  A  concise  way  to 
describe  such  wave  functions  is  that  they  are  the  unit  vectors  in  the  (free)  complex  vector  space 
on  the  set  (0,  1},  where  the  elements  of  (0,  1}  are  assumed  to  be  an  orthonormal  basis.  We  use 
Dirac  notation  to  denote  the  elements  of  this  basis  by  |0)  and  |1).  Thus  a  state  of  a  qubit  variable 
is  given  by  a  vector  a|0)  +  6|1)  where  a  and  h  are  complex  numbers  satisfying  \a\2  +  |6|2  =  1. 
Operationally,  when  measure  q  in  v  is  executed,  where  the  state  of  q  is  a|0)  +  b 1 1),  v  is  assigned 
0  with  probability  |a|2  and  is  assigned  1  with  probability  |6|2.  Suppose  that  in  a  particular  program 
P,  the  qubit  variables  are  q, ,  . . . ,  qn.  The  states  of  all  of  these  variables  in  the  aggregate  are 
the  unit  vectors  in  the  (free)  complex  vector  space  on  the  set  of  0, 1-valued  n-tuples,  which  is 
the  tensor  product  of  the  vector  spaces  corresponding  to  the  individual  qubits,  a  2n-dimensional 
complex  space.  The  latter  presents  a  complication:  take,  for  example,  two  qubits  p  and  q  jointly 
entangled  in  the  Bell  state  ^(|0)  <g)  |1)  +  |1)  <g)  |0)).  There  are  no  states  of  the  individual  qubits 
p  and  q  whose  tensor  product  is  the  Bell  state.  The  general  state  of  the  two  cubits  has  the  form 
«oo(|0)  <8>  |0»  +  a0i(|0)  <8>  |1))  +ai0(|l)  <8>  |0))  +an(|l)  <g)  |1)).  Should  a  measurement  of,  say  p  be 
performed  and  return  the  result,  say  1,  the  state  of  the  pair  jumps  to  (|ai0|2  +  |an|2)-1/2(aio(|l)  <8> 
|0))  +  an  ( 1 1)  ®  |  !)))•  Similar  jumps  happen  with  measurements  of  a  single  qubit  among  n  qubits. 
For  this  reason  the  state  of  all  of  the  qubit  variables  must  evolve  in  the  aggregate  during  the  course 
of  a  computation  and  explains  the  approach  we  take  below  to  the  qubit  components  of  the  heterotic 
dynamical  systems  we  employ. 

A  heterotic  dynamical  system  (greater  generality  is  possible  and  needed  for  some  applications 
than  we  give  here)  is  given  by  a  coupled  system  of  first  order  differential  equations 

{ ^  |  i  =  1,  ••.,  n 

where  each  Xi  is  a  function  from  the  reals  to  a  convergence  space  X,.  For  this  system  of  equations 
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to  make  sense,  the  derivatives  of  the  Xi  need  to  be  defined.  Furthermore,  to  capture  probabilistic 
nondeterministic  operations  such  as  occur  in  Abramsky’s  language,  we  must  have  available  among 
the  fi  “multi-valued  functions”  that  probabilistically  return  values  from  a  finite  set  of  possibilities 
[16]. 

The  sequential  nature  of  Abramsky’s  language  is  captured  in  a  heterotic  dynamical  system  by 
the  functions  fi  that  can  constrain  state  transitions  to  occur  only  at  specific  times  and  including 
an  auxiliary  variable  to  serve  as  a  program  counter.  For  assignment  commands  involving  integer 
arithmetic  expressions,  the  right-hand  side  of  the  corresponding  differential  equation  must  deter¬ 
mine  how  the  value  of  the  variable  is  to  change.  To  this  end,  the  set  of  integers  Z  is  equipped  with 
the  convergence  structure  corresponding  to  the  final  segment  topology  and  the  space  D(Jl,  Z )  of 
differentials  of  functions  from  the  reals  to  the  integers  is  chosen  to  consist  of  step  functions  of  a 
somewhat  special  kind:  for  each  pair  of  integers  m  and  n,  with  m  <  n,  the  step  functions  defined 
by  9t0  (t)  =  if  t  <  t o  then  n  else  m  and  hto  (t)  —  if  t  <  t0  then  m  else  n.  A  unitary  operation  U  on 
quantum  registers  (i.e.  the  qubits  taken  in  the  aggregate)  is  implemented  by  having  the  right-hand 
side  of  the  differential  equation  governing  the  register  produce  the  unitary  operator.  Regarding  pro¬ 
grams  in  Abramsky’s  language  as  heterotic  dynamical  systems  also  allows  to  couple  the  programs 
to  ordinary  real  and  complex  dynamical  systems  whose  states  evolve  in  continuous  time. 

In  classical  computation,  there  are  formal  verification  methods  based  on  specification  logic  [1]. 
Formulas  in  specification  logic  take  the  form 

mew 

where  formulas  $  and  T  are  expressed  in  an  underlying  logic  (usually  first-order  predicate  logic) 
and  where  C  is  a  command;  i.e.  a  program  or  fragment  of  a  program.  Formal  semantics  are  given 
by  the  notion  that  the  state  map  defined  by  C  maps  a  computation  state  about  which  $  is  true  to  a 
state  about  which  T  is  true;  i.e. 

[Wiwc  w 

where  [{C}]  is  the  map  on  computation  states  defined  by  the  command  C.  The  analogous  speci¬ 
fication  formulas  for  quantum  computation  can  be  taken  to  be 

{Trace  pP^  >p}U  {Trace  p  P^  >  q } 

where  A  is  an  observable  and  P£  =  xa(A),  the  characteristic  function  of  effective  Borel  set  A. 
(See  [21]  for  the  derivation  of  the  probability  that  the  result  of  a  measurement  of  observable  A 
with  be  in  A  if  a  system  is  in  a  quantum  state  given  by  positive  definite  trace  1  operator  p.)  U  is  a 
command  implementing  a  unitary  transformation 

An  effective  Borel  set  is  defined  as  follows.  Let  P(a,  . . . ,  xn)  be  a  computable  (1,  n)-ary 
relation,  where  the  X{  are  nonnegative  integers  and  cr  is  a  function  mapping  nonnegative  integers 
to  nonnegative  integers.  Computable  means  that  with  an  oracle  for  computing  a,  we  need  only 
use  the  oracle  to  compute  some  finite  sequence  a(0),  a(l),  . . . ,  a(k).  to  then  effectively  decide 
whether  or  not  P(a,  xi,  . . . ,  xn)  is  true.  Consider 

Vxi  3x2  ■  ■  ■  Qxn  P(a,  xi,  ,  xn) 

where  Q  is  either  V  or  3  in  the  alternating  quantifier  prefix  in  the  above  expression.  In  the  conven¬ 
tional  notation  widely  used  in  mathematical  logic,  the  sets  of  all  a  satisfying  formulas  such  as  the 
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above  are  the  (1,  0)-ary  arithmetic  relations  when  interpreted  on  the  natural  numbers  or  any  other 
effectively  presented  countably  infinite  set  such  as  the  integers.  These  sets  are  all  Borel  and  the 
set  of  their  characteristic  functions  are  the  Borel  sets  of  configurations  we  have  been  discussing. 
In  a  set-theoretic  sense  these  constitute  a  countable  subcollection  of  the  uncountably  many  distinct 
Borel  sets  of  configurations  we  have  been  implicitly  considering,  but  they  include  all  Borel  sets  of 
configurations  that  can  be  explicitly  defined  via  the  apparatus  of  formal  number  theory. 


5  CONCLUSIONS 

The  investigations  discussed  in  this  report  began  with  the  goal  of  developing  persistent  homology 
to  be  a  new,  additional  tool  for  software  validation.  After  an  initial  exploration  of  the  utility  of 
persistent  homology  for  detection  of  dynamic  geometric  structure  in  the  deformations  of  low  di¬ 
mensional  manifolds,  this  phase  of  the  project  concluded  that  applications  of  persistent  homology 
to  high  dimensional  complicated  dynamical  systems  is  premature  [5,  15].  The  project  then  fo¬ 
cused  on  developing  the  mathematically  rigorous  specification  and  formal  semantics  of  software 
involving  variables  with  compound  type  and  that  potentially  aggregates  more  than  one  model  of 
computation,  which  we  modeled  by  heterotic  dynamical  systems.  This  focus  on  rigorous  formal 
semantics  of  programs  for  complex  cyber  systems  was  judged  to  be  necessary  for  developing  for¬ 
mal  validation  methods.  The  capabilities  developed  during  this  investigation  have  provided  a  way 
forward  to  the  following  milestones: 

1 .  Development  of  a  specification  logic  for  verification  of  closely  coupled  classical  and  quan¬ 
tum  computational  processes. 

2.  Development  of  integration  on  convergence  spaces  for  rigorous  reasoning  about  the  compu¬ 
tation  of  the  trajectories  of  coupled  digital/analog  heterotic  systems. 

3.  Enhancement  of  formal  methods  for  synthesizing  programs  for  adiabatic  quantum  computa¬ 
tion. 
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7  LIST  OF  SYMBOLS,  ABBREVIATIONS  AND  ACRONYMS 


AFRL 

CONV 

CS 

PI 

PostDISCRETE 

PreTOP 

RDIGRAPH 

RITA 


Air  Force  Research  Laboratory 

The  Cartesian  closed  category  of  convergence  spaces 

Convergence  Structure 

Principle  Investigator 

The  category  of  post  discrete  convergence  spaces 
The  category  of  pretopological  spaces 
The  category  of  digraphs  with  self-loops 
Emerging  Computing  Architectures 
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